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The bacterial potassium channel KcsA is gated open by the binding of protons to amino acids on the intracellular 
side of the channel. We have identified, via channel mutagenesis and x-ray crystallography, two pH-sensing amino 
acids and a set of nearby residues involved in molecular interactions that influence gating. We found that the mini- 
mal mutation of one histidine (H25) and one glutamate (E118) near the cytoplasmic gate completely abolished 
pH-dependent gating. Mutation of nearby residues either alone or in pairs altered the channel's response to pH. 
In addition, mutations of certain pairs of residues dramatically increased the energy barriers between the closed 
and open states. We proposed a Monod-Wyman-Changeux model for proton binding and pH-dependent gating 
in KcsA, where H25 is a "strong" sensor displaying a large shift in pK a between closed and open states, and El 18 is 
a "weak" pH sensor. Modifying model parameters that are involved in either the intrinsic gating equilibrium or the 
pK a values of the pH-sensing residues was sufficient to capture the effects of all mutations. 
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INTRODUCTION 



Ligand-gated ion channels convert the energy of ligand 
binding into the gating of ion-conducting pores, allow- 
ing the passage of current across biological membranes. 
In this way, electrical activity under the explicit control 
of signaling molecules underlies many physiological pro- 
cesses such as postsynaptic neurotransmission in nerve 
and muscle, pain sensation, and sensory transduction 
(Hille, 2001). 

The proton-gated prokaryotic K + channel KcsA 
(Schrempf etal., 1995; Cuello etal., 1998; Heginbotham 
et al., 1999) was the first ion channel modeled with 
high resolution x-ray crystallography (Doyle et al., 1998; 
Zhou et al., 2001). The pore appeared closed, as the 
inner TM2 helices formed a tight hydrophobic constric- 
tion called the "bundle crossing." Spectroscopic studies 
confirmed that the conformation of the bundle-cross- 
ing gate was pH dependent (Perozo et al., 1998; 
Liu et al., 2001). Structural studies of MthK (Jiang 
et al., 2002a), a Ca 2+ -gated K + channel (Zadek and 
Nimigean, 2006; Li et al., 2007), provided a model for 
the open conformation (Jiang et al., 2002b), later found 
to be consistent with the presumed open-state struc- 
ture of a mutant KcsA pore (Cuello et al., 2010b). These 
open versus closed structural snapshots represent the 
conformational states of KcsA that equilibrate in a pro- 
ton-dependent manner. 
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In a previous report, we made several mutations near 
the bundle-crossing gate of KcsA that dramatically al- 
tered pH-dependent gating (Thompson et al., 2008). 
In particular, we confirmed a contribution from a key 
histidine (H25) behind the bundle crossing that had 
been proposed as the sole pH sensor in a nuclear mag- 
netic resonance (NMR) study (Takeuchi et al., 2007). 
Despite a significant effort by us and others to modify 
or ablate pH sensing using mutagenesis (Thompson 
et al., 2008; Cuello et al., 2010a,c; Hirano et al., 2011), 
the energetic mechanism underlying proton-dependent 
gating remained obscure partly because of ambiguity 
of interpretation after channel modification. Muta- 
tions may alter the energy of ligand binding, the ligand- 
independent intrinsic energy of gating, or both, and 
these outcomes are often not distinguishable by analy- 
sis of channel activity alone (Colquhoun, 1998). Further- 
more, there is a high density of arginines and ionizable 
residues at this location, and mutations of these resi- 
dues may change proton binding (pK a ) to nearby pH- 
sensing sites by alteration of local electrostatics or 
solvent accessibility, making the determination of pH- 
sensor contributions difficult. 

In this study, we investigated the specific role of each 
ionizable residue located at the bundle crossing of KcsA 
using channel recordings in lipid bilayers and x-ray 
crystallography of two bundle-crossing mutant chan- 
nels. We propose that proton binding to H25 and El 18 
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is energetically linked to proton-dependent gating in 
KcsA, and thatR121, R122, and E120 play a supportive 
role. We propose that specific intersubunit and intrasu- 
bunit interactions alter the intrinsic equilibrium between 
closed and open states, whereas proton binding gates 
the channel according to a Monod-Wyman-Changeux 
(MWC) kinetic scheme. In our model, H25 is a "strong" 
pH sensor, characterized by a large pK a shift between 
closed and open conformations, enabling this residue 
to change the gating equilibrium to a large degree. On 
the other hand, El 18 is modeled as a "weak" pH sen- 
sor, suggesting a smaller energetic contribution from 
this residue. 



MATERIALS AND METHODS 

Protein expression, purification, and reconstitution 
in lipid bilayers 

KcsA mutations were introduced with the QuikChange site- 
directed DNA mutagenesis kit (Agilent Technologies) and verified 
with sequencing. All pH-sensor mutations were made on the back- 
ground of the E71A mutant, shown to remove selectivity filter 
inactivation (Cordero-Morales et al., 2006). KcsA protein was ex- 
pressed and purified as described previously (Heginbotham et al., 
1999; Thompson et al., 2008). All reagents were from Sigma- 
Aldrich unless otherwise noted. In brief, either N-terminal 6His- 
tagged KcsA pASK90 or C-terminal 6His-tagged KcsA pQE60 
constructs were transformed into JM83 (ATCC) or BL21 (DE3; In- 
vitrogen) Escherichia coli cells, respectively, and grown at 37°C in 
Terrific Broth or Luria-Bertani media to OD 60 o = 1- For pASK90, 
protein expression was induced for 90 min with 0.2 mg/ml anhy- 
drotetracycline (Acros Organics), and for pQE60, protein expres- 
sion was induced for 3 h with 500 pM IPTG. Because of poor 
expression at 37°C, H25R/E118A expression was induced at 20°C 
overnight. Cells were harvested via centrifugation at 5,000 gfor 
15 min at 4°C, resuspended in 100 mM KC1 and 50 mM Tris, pH 7.5, 
and disrupted with probe sonication (Thermo Fisher Scientific). 
Membranes were extracted for 2 h at room temperature with 
50 mM n-decyl maltoside (DM; Affymetrix) followed by centrifu- 
gation at 35,000 gfor 45 min. Protein was purified from the solu- 
ble fraction in buffer (100 mM KC1, 20 mM Tris, and 5 mM DM, 
pH 7.5) with either a Ni 2+ -affinity column (EMD Millipore) or 
a HiTrap Co 2+ -affinity column (GE Healthcare), eluted with the 
addition of 300 mM imidazole, and finished with gel filtration 
(Superdex 200; GE Healthcare). 

Purified channel was reconstituted into liposomes at protein to 
lipid ratios of 0.1-10 pg protein per milligram of lipid (3:1 POPE/ 
POPG; Avanti Polar Lipids, Inc.). Detergent was removed using a 
G50 fine (GE Healthcare) gel-filtration column in buffer contain- 
ing 400 mM KC1, 5 mM NMG, and 20 mM Tris, pH 7.5. Liposome 
aliquots were flash-frozen in liquid nitrogen and stored at — 80°C. 

Single-channel recording and analysis 

KcsA pH-sensor mutants were recorded and analyzed as described 
previously (Thompson et al., 2008). In brief, channels were re- 
corded in a horizontal planar lipid bilayer system consisting 
of two chambers separated by a plastic partition (transparency 
slide; IKON) containing a 50-100-pm hole, across which bilayers 
were formed using a mixture of POPE (7.5 mg/ml) and POPG 
(2.5 mg/ml) dissolved in rc-decane. Recording solutions were: 
70 mM KC1, 30 mM KOH, 10 mM MOPS, 10 mM succinate, and 
10 mM Tris, pH adjusted to 4-9. Proteoliposomes were applied to 
the cis (upper) chamber, which contained recording solution at 



pH 7. Initially, the trans (lower) chamber contained recording 
solution at pH 4, such that only channels with pH sensor facing 
trans were activated. Only channels that displayed typical KcsA 
behavior in these conditions were studied. The recording elec- 
trode was placed in the cis chamber, and the ground electrode was 
placed in the trans chamber. 

Proton dose-response curves were measured by perfusing the 
trans chamber with at least 5 ml of recording solutions with vary- 
ing pH and measuring channel activity. Bilayers were voltage- 
clamped using an amplifier (Axopatch 200A, 200B; Molecular 
Devices) or a bilayer clamp (BG535; Warner Instruments). Single- 
channel currents were filtered at 2 kHz using a four-pole Bessel 
filter, digitized at 25 kHz using a Digidata 1440A or 1320 (Molecular 
Devices), and recorded in Clampex 10 or Axoscope 10 software 
(Molecular Devices). Channel orientation was confirmed by re- 
cording the typical conductance and open-channel noise proper- 
ties of KcsA E71A. Single-channel current amplitudes and open 
probabilities (Po) were analyzed in Clampfit 10 software (Molecu- 
lar Devices). The number of channels in the bilayer during an 
experiment was verified by perfusing back to pH 4.0 during the 
experiment. Bilayers where this was not achieved (because of 
bilayer rupture or gain/loss of channels) were not included in 
the final analysis. 

Po versus pH data (not normalized) were plotted in Origin (Mi- 
crocal) and fit with the model-independent Hill equation (Eq. 1 ) : 
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where P 0 max is the maximal open probability, EC 50 is the proton 
concentration of half-activation (also represented as pHi /2 ) , [If] 
is the proton concentration (also represented as pH), and n H is 
the Hill coefficient, which describes the steepness of the dose-re- 
sponse curve. 

The Po versus pH data were also fit to either a single proton- 
binding site MWC model (Fig. 4 B and Eq. 2) : 
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where A is given by: 
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or a two proton-binding site MWC model (Fig. 4 C and Eq. 3) 
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In these models, proton binding was assumed to be equivalent 
and independent to all KcsA subunits. L a is the intrinsic gating 
equilibrium between open and closed in the absence of protons, 
[ft] is the proton concentration (also represented as pH), and 
j^open an£ j j^dosed are ac y dissociation constants for a proton 
sensor in the open and closed states, respectively (also repre- 
sented as pK° p,!n and pK a dosed ) . Eq. 3 contains two proton sensors 
that are labeled with subscripts 1 and 2. In general, the absence of 
direct experimental constraint on the values for L Q and the very 
steep [iL] dependence for many of the datasets invalidated de- 
velopment of uniquely determined model parameters for individ- 
ual datasets. We made the simplifying assumption that mutation 
of pH-sensor site 1 (identified as H25) did not change the pK a 
values for pH-sensor site 2 (identified as El 18), and vice versa. In 
addition, we also assumed that mutation of E120 did not perturb 
nearby pH-sensor pK a values. Collectively, these assumptions al- 
lowed us to group several datasets together and provided signifi- 
cant constraint on the pK a values that could describe the data. We 
globally fit all mutations of H25 that retained the E 1 1 8 pH sensor. 
Mutants H25R and H25R/E120A were fit simultaneously to a 
single pH-sensor model (El 18) with shared values for pK^^" and 
/?X„ closed but with two separate values for L>. The fit was unique 
and well-determined but with large uncertainty in the L> values 
^closed, Eiis = g g ± 0 1; p K °p™. E»8 = 6.2 ± 0.1, L 0 H23R = 7 x 1(T 4 ± 
4 x lfT 4 , and l o H25R/e120A = 0.08 ± 0.02). Next, we attempted to 
globally fit all mutations of El 18 that retained the H25 pH sensor. 
Mutants E118A, E118Q/E120Q, and E118A/E120Awere simulta- 
neously fit to a single pH-sensor model with shared values for 
^closed and pK°p™ but t h ree separate L, values. In this case, 
the pK a values were not uniquely determined but varied over a 
range of M).5 pH units or less as parameter initialization was sys- 
tematically varied. Therefore, the approximate pK a values were 
fairly well constrained by the requirement to simultaneously fit 
both the E118A and E118A/E120A dose-response curves that 
were shifted by ~1.5 pH units. We found a global fit that overall 
minimized the fitting error but with large uncertainty remaining 
in L, values {pK***^ H25 = 4.8 ± 0.6, pKf?™- H25 = 7.6 ± 0.7, L, E118A = 
HT 8 ±8x 10- 8 ,Z c E118 ^ E120 «=0.8x KT 8 ±5.8x 10- 8 ,andL o E118A/E120A = 
2 x 10~ 4 ± 9 x 10~ 4 ). The pK a values for H25 were subsequently 
constrained to these best estimates, and the Lo values for all mu- 
tants were fitted using those constraints (these values are reported 
in Table 1 ) . A similar procedure was followed for the unknown 
pH sensor in the triple mutants H25R/E118A/E120A (Fig. 1), 
H25R/E118Q/E120Q, and H25A/E118A/E120A (Fig. S2). Only 
steady-state Po data were fit to models (Figs. 1-3). No kinetic analy- 
sis of single-channel data was performed. All data fits are summa- 
rized in Table 1 . 

The change in the energy of gating caused by a single proton 
binding to a pH sensor is given by Eq. 4: 
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The gas constant, R, times the temperature T in Kelvin is RT = 
0.59 kcal/mole at 25°C, and the coupling constant, 7, is given by, 
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For the two proton-binding site model, (3 is used for the second 
site-coupling constant. 

The pH-independent equilibrium constant for the double 
mutant H25R/E118A was predicted from our models using the in- 
trinsic equilibrium constants for the control channel and the single 
mutants H25R and E118A by, 
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Crystallography 

As described previously (Zhou et al., 2001), the C-terminal do- 
main of the KcsA El 18A or the KcsA R122A mutant was removed 
using chymotrypsin before gel-filtration chromatography to aid 
in crystallization. The purified channel was then combined with a 
KcsA antibody F ab fragment and run over a second gel-filtration 
column (Zhou et al., 2001). The KcsA-F ab complex was concen- 
trated to 8-10 mg/ml. Protein was crystallized in hanging drop 
vapor diffusion trays by combining 2 pi of protein solution with 
2 pi of a well solution consisting of 19-22% (wt/vol) PEG 400, 
25 mM Mg acetate, and 100 mM MES, pH 6.5. After crystals 
formed, in about 2 wk, the PEG 400 concentration in the well 
solution was adjusted to 40% for cryoprotection. After 24 h, the 
crystals were flash-frozen in liquid nitrogen. 

Diffraction data were collected at 93 K using the X25 beamline 
at the National Synchrotron Light Source at Brookhaven National 
Laboratory. Reflections were indexed, integrated, and scaled in 
space group 14 using HKL2000 (HKL Research). The data were 
phased using molecular replacement with Protein Data Bank ac- 
cession number 1K4C as the starting model. Phasing, refinement, 
and structure quality assessment were performed with the Phenix 
suite of programs (Adams et al., 2010). Model building was imple- 
mented with Coot (Emsley et al., 2010). Models and structure 
factors for the mutant KcsA structures E118A and R122A are de- 
posited in the Protein Data Bank under accession numbers 4LCU 
and 4LBE, respectively. 

Online supplemental material 

Fig. SI is provided to illustrate the modest changes in KcsA 
gating kinetics for the individual mutations E118A, R121A, 
and R122A for comparison, with the results presented in Fig. 2. 
Fig. S2 summarizes evidence for the existence of unidentified 
pH sensors in two mutants with the H25 and El 18 sensors re- 
moved. Fig. S3 shows the pH dependence of the H25A mutant, 
demonstrates the dramatic channel to channel variability in this 
mutant, and further strengthens the argument for the necessity to 
model at least two pH-sensor residues. The online supplemental 
material is available at http://www.jgp.org/cgi/content/full/ 
jgp. 20131 1057/DC1. 



RESULTS 

Double mutation of histidine 25 and glutamate 
118 ablated pH-dependent gating 

The KcsA proton-dependent gate was identified as the 
cytoplasmic constriction formed by the bundle crossing 
of inner transmembrane helices (Perozo et al., 1998; 
Gross et al., 1999; Kelly and Gross, 2003). Several titrat- 
able residues are proximal to this channel gate (Fig. 1 A) , 
indicating a possible role in pH sensing. In our previous 
report, we combined the mutation of two such glutamates, 
E118A/E120A, with the H25R mutant and established 
a near-pH-insensitive channel (Fig. 1 B) , suggesting 
that the major components of pH sensing were gone 
(Thompson et al., 2008). H25 was implicated as a pH 
sensor because of a significant reduction in the steep- 
ness of the proton dose-response described by the Hill 
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coefficient. El 18/E120 were implicated as pH sensors 
because of a large shift relative to control in pH de- 
pendence from ^5 to ^6.5 pH (Fig. 1 B) (Thompson 
etal., 2008). 

We investigated whether or not both glutamates con- 
tributed to pH-dependent gating as important proton- 
binding sites (Fig. 1, C and D). As in our previous report 
(Thompson et al., 2008), all channel mutations were 
made on the background of E71A, a mutation that re- 
moves selectivity filter inactivation and reveals full 
activation of the pH-dependent gate at low pH values 
(Cordero-Morales et al., 2006) . In this paper, KcsA E71A is 
referred to as "KcsA" or "control channel." H25R/E120A 
remained pH sensitive with a shifted proton dependence 
relative to the H25R mutant. H25R/E118A on the other 



hand, exhibited no pH dependence and full activation 
in all pH conditions (Fig. 1, C and D). Example single- 
channel recordings of E118A, H25R, and H25R/E118A 
are shown to illustrate the combined contribution from 
these two residues during pH-dependent gating (Fig. 1 E) . 
Importantly, H25R/E118A resulted in very little chan- 
nel to channel variability in contrast to the previously re- 
ported H25R/E118A/E120A mutant (Thompson et al, 
2008). These results suggested that E120 may not directly 
contribute to pH-dependent gating as a critical proton- 
binding site but still contributes to the gating equilibrium 
as seen in the double mutants E118A/E120A and H25R/ 
E120A (Fig. 1, B and C). On the other hand, proton 
binding to both H25 and El 18 may account for the full 
pH-dependent activation of KcsA, assuming that the 
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Figure 1 . Protonation of H25 
and El 18 may be sufficient to 
fully activate KcsA. (A) A high 
density of titratable and charged 
residues surrounds the cytoplas- 
mic gate of KcsA, viewed from 
the side (left) and from the cyto- 
plasm (right). H25, El 18, E120, 
R121, and R122 are highlighted. 
(B) Summary of pH-dependent 
open probability (Po) for critical 
mutants reported in Thompson 
et al. (2008). Lines are fits with 
MWC models (Materials and 
methods; Fig. 4 and Table 1), 
where H25 is proton-binding 
site 1 and the El 18 is site 2. The 
E71A control channel (black 
squares) was fit to a two pH-sen- 
sor model (black dashed line). 
H25R/E118A/E120A retained a 
minor residual pH dependence 
(open triangles, black line) , mod- 
eled as a single proton-binding 
site. (C) E120 is likely not a pH 
sensor, and E120A was mod- 
eled similar to the control chan- 
nel with two pH-sensor residues 
(circles, orange). H25R/E120A 
retained the El 18 pH sensor 
(squares, blue line). (D) E118A 
was modeled with H25 pH sen- 
sor (circles, green line). H25R/ 
E118A (squares, black line) was 
fitted to a pH-independent gat- 
ing equilibrium (L„ = 51). Data 
are the mean ± SEM from three 
to seven bilayers, and the model 
for wild type is shown as a black 
dashed line in C and D. (E) Ex- 
ample single-channel recordings 
for E 1 1 8A ( top) , H25R (middle ) , 
and H25R/E118A (bottom) at 
three pH values illustrating the 
pH-dependent Po. 10 s of data at 
100 mV were filtered offline to 
500 Hz for display. 
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mutations (H25R/E118A) did not alter the intrinsic 
gating of the channel through steric or structural changes 
that do not mimic protonation in wild type. 

Individual mutations of the glutamates showed that 
although E120A was nearly indistinguishable from the 
control channel (Fig. 1 C), the El 18A mutation (Fig. 1 D) 
also resulted in only a very small shift in pH dependence, 
despite the dramatic effect of this mutation when com- 
bined with H25R. Considering the small effect observed 
when only El 18 is mutated compared with the consider- 
able effect of the H25R mutation, we propose that pro- 
tonation of both El 18 and H25 is necessary to gate the 
channel, and that H25 contributes a considerably greater 
amount of energy during the gating process. 

Two arginine residues influence pH dependence 
and gating kinetics 

We next investigated the role of residues near the pro- 
posed El 18 and H25 pH sensors. Two arginine residues, 



R121 and R122, appear to interact with El 18 in the 
closed-state KcsA crystal structure (Zhou et al., 2001; 
Thompson et al., 2008). We explored the effect these 
residues have on pH-dependent gating by mutating them 
individually to alanine. R121A introduced a significant 
shift in pH dependence (Fig. 2 A, blue squares) , with no 
change in steepness (solid blue line) relative to control 
(black dashed line). The R122A mutation (Fig. 2 B, 
green squares) resulted in a relatively modest shift in 
pH dependence as well as a small decrease in slope 
(n H = 3.4; Table 1). 

We next removed each individual arginine on the 
background of the E118A mutation (Fig. 2, A and B). 
These mutant channels displayed gating kinetics that 
differed remarkably from the control channel. Fig. 2 C 
illustrates the kinetic properties observed for these 
channels at intermediate pH values where the channels 
are near half-activation. Unlike the control channel, 
which opened and closed on the millisecond timescale, 
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may contribute to the in- 
trinsic gating equilibrium as 
well as to the pK a values of 
nearby proton-binding sites. 
(A) R121A Po versus pH re- 
lationship (blue squares) was 
modeled either by changing 
relative to control only the 
intrinsic gating (solid blue 
line; model fit 1 in Table 1) 
or by also changing pK a val- 
ues (dashed blue line; model 
fit 2 in Table 1). The E118A/ 
R121A pH dependence (red 
circles) was not accurately 
25 s ~1 min. -1.2 min. established, and the red dot- 

ted line has no theoretical 
meaning. (B) R122A Po ver- 
sus pH relationship (green 

_5 m j n -12 min squares) was modeled either 

by changing relative to con- 

Q trol only the intrinsic gating 

(solid green line; model fit 1 
in Table 1 ) or by also chang- 
ing pK a values (dashed green 
line; model fit 2 in Table 1). 
The E118A/R122A pH de- 
pendence (orange circles) was 
not accurately established, 
and the orange line has no 
theoretical meaning. Data in 
A and B are the mean ± SEM 
from three to eight bilayers, 
and the model for control is 
shown as a black dashed line. 

MWC model fit parameters used in A and B are given in Table 1. (C) Example recordings of E71A control, E118A/R122A, and E118A/ 
R121A at 100 mV and filtered offline to 50 Hz (or 20 Hz for E118A/R121A) for display. (D) The intracellular end of the TM2 helix (gray 
ribbon) of the El 18A KcsA mutant aligned with the same region from KcsA wild type (Protein Data Bank accession no. 1K4C) . The resi- 
dues of interest are displayed as red sticks (E118A mutant) or with carbons in yellow, nitrogens in blue, and oxygens in red (wild type). 
The electron density map for the El 18A KcsA mutant is shown contoured at 1.2 a (Table 2) . R121 was more disordered compared with 
wild type, and R122 assumed a new conformation. 
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E118A/R122A displayed minute-long open bursts and 
closed durations. Even more dramatically, the E118A/ 
R121A mutant displayed very long stretches of high 
activity (^10 min) interrupted by equally long closings. 
Brief closures within a burst, called channel flickers, are 
thought to arise from residual selectivity filter gating 
events (Piskorowski and Aldrich, 2006) and are likely 
not related to movement of the pH-dependent gate. 
Although both mutants were clearly pH dependent, the 
dose-response curves were difficult to establish accu- 
rately because of the uncertain Po values measured at 
intermediate pH values (Fig. 2, A and B, red and orange 



circles) . Regardless, these data have interesting implica- 
tions. First, they show that residues within the pH-sensing 
region are able to dramatically increase the energy bar- 
riers between the open and closed states, resulting in 
drastically slower gating. These intermediate states dur- 
ing gating likely relate to the intersubunit cooperativity 
and may also explain why an increase in subconduc- 
tance gating is also observed in these mutant channels 
(Fig. 2 C). Second, the individual mutations E118A and 
R121A both slow channel gating modestly within the 
time scale of seconds, whereas R122A has very little ef- 
fect on gating kinetics (Fig. SI). The drastic change in 
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(0.1) 


(0.1) 


(4 x 10"") 




0.3 


1 9 
1.3 






J.yf 


K 9 


U.Uo 




(0.04) 


(0.1) 






(0.1) 


(0.1) 


(0.02) 


El 18A 


O.O-i 

(0.03) 


-i.J 

(2.8) 


4 8 


7 6 






1 1 x 10 s 
(0.1 x 10" 8 ) 


V 1 1 SO /F 1 90O b 


j.j 

(0.004) 


4 4 
(1-2) 


4 8 


7 6 






8 1 x 10~ 9 
(0.5 x 10"') 


E118A/E120A b 


6.6 

(0.003) 


5.8 
(0.1) 


4.8 


7.6 






1.5 x 10" 4 
(0.4 x 10" 4 ) 


H25R/E118A 














52 
(11) 


E120A 


5.3 

(0.006) 


3.4 
(0.1) 


4.8 


7.6 


5.0 


6.2 


2.5 x 10" 12 
(1.2 x 10" 12 ) 


H25R/E118A/E120A b - d 










4.7 


6.5 


3.3 
(0.7) 


H25R/E118Q/E120Q M 


6.3 

(0.04) 


2.5 
(0.5) 






4.7 


6.5 


0.021 
(0.004) 


H25A/E1 18A/E120A b ' d 


6.3 
(0.01) 


2.6 
(0.1) 






4.7 


6.5 


0.014 
(0.002) 


R121A (model fit 1) 


5.8 
(0.01) 


4.0 

(0.2) 


4.8 


7.6 


5.0 


6.2 


1.3 x 10"' 
(0.3 x 10" 9 ) 


R121A (model fit 2) 






6.6 


7.6 


5.0 


7.7 


1.0 x 10" 11 
(0.1 x 10" 11 ) 


R122A (model fit 1) 


5.6 
(0.01) 


6.9 
(0.7) 


4.8 


7.6 


5.0 


6.2 


3.8 x 10" 11 
(0.6 x 10" 11 ) 


R122A (model fit 2) 






6 


7.6 


5.0 


6.2 


4 x 10" 9 

(io- 9 ) 


E118A/R121A (not modeled) 


6.0 

(0.2) 


0.6 

(0.2) 












E118A/R122A (not modeled) 


6.6 

(0.09) 


1 .3 

(0.4) 












E120A/R122A 


5.7 

(0.02) 


2.1 

(0.2) 


6 


6.5 


4.5 


5.5 


4.2 x 10~ 3 
(0.3 x 10~ 3 ) 



Values in parentheses are the standard error of fit. 
a All mutants on the background of E71A. 
b Data and Hill fits from Thompson et al. (2008). 

c Not applicable; no error given means the parameter was constrained {see Materials and methods) . 
d Modeled as an unidentified residual pH sensor. 
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the gating kinetics observed for both double mutants 
(E118A/R121Aand E118A/R122A) indicates energetic 
coupling interactions between El 18 and both R121 and 
R122 during gating transitions. 

We further investigated the structural relationship 
between E118 and the two arginines, R121 and R122, 
by solving the crystal structure of the E118A mutant in 
the closed state (Fig. 2 D and Table 2). Supporting our 
conclusion that El 18 interacts with both arginines, the 
structure of E118A shows deviations from the wild-type 
structure for both R121 and R122. R122 adopted a 
different conformation, whereas R121 became disor- 
dered in the mutant structure (Fig. 2 D, residues shown 
in red). The mutation also resulted in the disruption 
of the potential intersubunit R122-E120 interaction 
(Fig. 1 A) . The interactions these residues make in the 
open state and during gating intermediate states re- 
main unknown. 

An intersubunit salt bridge does not contribute 
significantly to pH sensing 

We proposed previously that an intersubunit salt bridge 
between R122 and El 20 observed in the closed-state 
KcsA crystal structure contributes to pH-dependent gating 



(Thompson et al., 2008). However, as shown above, the 
E120A mutation alone does not alter pH-dependent 
gating (Fig. 1 C). Furthermore, E120A in addition to 
the R122A did not further perturb pH-dependent gat- 
ing (Fig. 3 A). A crystal structure of the R122A mutant 
showed almost no change in the structure of the bundle 
crossing and proposed pH-sensor residues (Fig. 3 B 
and Table 2). It is possible that R122 can make several 
molecular interactions that modestly alter the intrinsic 
closed-open equilibrium; however, these results fur- 
ther argue that El 20 is not a pH-sensor residue and that 
proton binding to El 18 and H25 alone drives gating. 

DISCUSSION 

Several previous studies have used mutagenesis targeted 
to modify or remove KcsA pH-dependent gating to es- 
tablish the proton-binding sites and qualitatively de- 
scribe the mechanism (Takeuchi etal., 2007; Thompson 
et al., 2008; Cuello et al., 2010c; Hirano et al., 2011). 
Hirano et al. (2011) proposed that the C-terminal do- 
main contains strong pH sensors at E146 and D149. 
However, truncation of the bottom half of the C-termi- 
nal domain (A140) (Cortes et al., 2001) as well as its 



TABLE 2 

Crystallographic data collection and refinement statistics 
KcsAE118A 



KcsAR122A 



Data collection 

Space group 
Cell dimensions 

a, b, c (A) 

<*. P,7(°) 
Resolution (A) 

I/o-I 

Completeness (%) 
Redundancy 



14 

155.95, 155.95, 75.95 
90.00, 90.00, 90.00 
50.0-2.75 
(2.85-2.75) a 
10.8 (93.2) 
18.5 (2.4) 
99.7 (99.3) 
7.4 (7.1) 



14 

155.94, 155.94, 75.81 
90.00, 90.00, 90.00 
50.0-2.75 
(2.85-2.75) 
8.3 (65.5) 

20.8 (2.2) 

99.9 (99.2) 
7.0 (6.2) 



Refinement 

Resolution (A) 
No. reflections 

Rwork/Rfree 

No. atoms 

Protein 

Ligand/Ion 

Water 

B factors 

Protein 

Ligand/Ion 

Water 

Root mean square deviations 

Bond lengths (A) 
Bond angles (°) 

a Values in parentheses are for the highest-resolution shell. 



49.31-2.75 

23,757 
0.183/0.216 

4,057 
37 
56 

83.6 
79.2 
57.4 

0.008 
1.481 



42.87-2.75 

23,755 
0.169/0.202 

4,081 
44 
SO 

72.9 
85.0 
56.7 

0.012 
1.819 
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replacement with a cytoplasmic domain from a differ- 
ent ion channel (unpublished data) did not change 
KcsA pH dependence, suggesting that this region has 
little effect on KcsA gating. An NMR study (Takeuchi 
etal., 2007) identified the importance of H25; however, 
no single mutation of H25 removed pH dependence in 
electrophysiological assays (Fig. 1 and Thompson et al., 
2008) or in radioactive flux assays (Cuello et al., 2010c) . 
We also examined H25A (the substitution used in the 
NMR study) to compare the results with H25R. Unfor- 
tunately, this mutant results in a high degree of channel 
to channel variability, but importandy, pH dependence 
remains (Fig. S3). Furthermore, single point mutations 
at other candidate residues (H20, El 18, E120, or H124) 
did not ablate pH dependence (Thompson etal., 2008; 
Cuello et al., 2010a). Therefore, functional analysis 
pointed to the likely contribution from at least two pro- 
ton-binding sites. 

A proposed model for proton binding and gating 
in the KcsA channel 

The proton-binding and gating relationship in KcsA ap- 
pears to be a highly cooperative process (n H = 4.4 for KcsA 
control; Table 1 ) . As in other ligand-gated channels, this 
property can arise naturally from subunit cooperativity 
during gating. We propose an MWC model for KcsA pH 
sensing and the interpretation of our mutagenesis re- 
sults. The fundamental assumptions underlying the MWC 
model are that the protein exists in two global conforma- 
tions, open (O) and closed (C), to which the ligand binds 
with different affinities (Monod et al., 1965; Changeux, 
2012). X-ray crystallography has been used previously to 
determine both closed and open conformations of KcsA, 
shown in Fig. 4 A (Zhou etal., 2001; Cuello etal., 2010b). 

The simplest pH-sensor MWC model with a single pro- 
ton-binding site per KcsA subunit yields a 10-state model 
with 5 closed and 5 open states, one each for the differently 



liganded states (Fig. 4 B). On the other hand, two inde- 
pendent proton-binding sites per KcsA subunit (as 
proposed above: H25 and El 18) can lead to a larger 
MWC model, with protonation occurring along two in- 
dependent directions, one for each proton-binding site 
(Fig. 4 C). For each differently liganded state in the 
single proton-sensor model, there are 5 differently li- 
ganded states for the secondary proton sensor, leading to 
a total of 25 closed and 25 open states. 

We thus formalized the contributions from the two 
pH-sensor residues (H25 and El 18) using a two pro- 
ton-binding site MWC model that emphasizes subunit 
cooperativity and state-dependent pK a values that ener- 
getically drive gating. All other contributions to the 
protein chemical energy combine to set the intrinsic 
equilibrium between open and closed conformations, 
called L Q . In this view, descriptions of interacting amino 
acids and electrostatic repulsion, variously emphasized in 
the literature (Miloshevsky and Jordan, 2007; Thompson 
et al., 2008; Cuello et al., 2010c), are connected to a spe- 
cific and testable chemical property, state-dependent 
pK a values. We used these models to describe the open 
probabilities (Po) versus pH (Eqs. 2 and 3; Materials 
and methods) of the KcsA control channel (containing 
the two proposed pH sensors, El 18 and H25) , the H25R 
mutant (gated by the single El 18 sensor), and the E118A 
mutant (gated by the single H25 sensor) . The individual 
datasets do not provide enough information to constrain 
the MWC model parameters pardy because the Po ver- 
sus pH rises steeply from 0 to 1 . Thus, to constrain the 
pK a values, we assumed that mutation of the H25 sensor 
did not perturb the pK a values of the El 18 sensor and 
vice versa. These assumptions allowed us to group our 
datasets into global fits that added a significant constraint 
on the pH-sensor pK a values (see Materials and methods) . 
Model fits are plotted in Figs. 1-3, and parameters are 
summarized in Table 1 . 




El 20 



R122 (wt) 



Figure 3. An intersubunit interaction be- 
tween R122 and El 20 is not essential for pH- 
dependent gating in wild type. (A) The pH 
dependence of E120A/R122A (red circles) 
does not significantly differ from the R122A 
mutant (blue squares). Two pH-sensor resi- 
dues are modeled for R122A (blue line) and 
E120A/R122A (red line), as in the control 
channel (black dashed line) . MWC model fit 
parameters are given in Table 1. (B) Align- 
ment of the intracellular end of the TM2 
helix (gray ribbon) of the wild-type KcsA crys- 
tal structure (Protein Data Bank accession no. 
1K4C; residues are modeled as sticks with car- 
bon yellow, nitrogen in blue, and oxygen in 
red) with a model of our R122A KcsA mutant 
structure (red sticks). The electron density 
map for E118A contoured at 1.2 cr (Table 2) 
showed that the mutation did not significantly 
alter the location of El 20 or the other nearby 
residues El 18, Q119, and R121. 
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Because mutation of H25 resulted in a significant re- 
duction of proton sensitivity (decrease in steepness of 
proton dose-response) and mutation of the El 18 did 
not, we modeled H25 as a strong proton sensor with a 
large pK a shift between closed (pK a closed = 4.8) and open 
(pK a °v en = 7.6). The single-site MWC model with these 
values approximated well the pH dose-response of the 
El 18A KcsA mutant that has just the H25 proton sensor 
remaining (Fig. 1 D). The glutamate proton sensor was 
described as weak, with a more modest pK a state depen- 
dence (pK a clos<!d = 5.0 and pK°^ n = 6.2). The single-site 
MWC model with these values also approximated well 
the pH dose-response of the H25R KcsA mutant that 
has just the El 18 proton sensor remaining (Fig. 1 B). 
For the models to fit the dose-response curves of the 
mutant channels, we assumed that mutation of a pH- 
sensor site did not alter the pK a values of the remaining 
pH sensor, but rather changed the intrinsic gating equi- 
librium constant, L a (Fig. 4 B and Table 1; see Materials 
and methods) . This made sense because the mutant 



channels were expected to have increased L Q values 
compared with the control channel, as the pH-sensor 
sites were mutated in a way to mimic constitutive pro- 
tonation, hence favoring channel opening. These mod- 
els, proposed as a conceptual starting point for further 
investigation, lack sufficient constraints to define the 
precise energetic contributions from the two proton 
sensors, El 18 and H25. However, we notably captured 
the salient feature in our data that mutation of H25 
always significantly decreased the slope of the pH dose- 
response, suggesting that the individual energetic con- 
tribution from H25 is significantly greater than the 
contribution from other proposed pH-sensor sites. 

In these models, the pH sensors H25 and El 18 have 
state-dependent pK a values, such that the affinities for 
proton binding are higher in the open state (higher pK a 
values) . This makes structural sense for H25, which ap- 
pears buried in the closed-state crystal structure (Zhou 
et al., 2001) and near three arginine residues (R117, 
R121, and R122), likely hindering proton binding and 
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Figure 4. The MWC model links proton 
binding to the equilibrium between closed 
and open KcsA conformations. (A) Crystal 
structures of closed KcsA (Protein Data Bank 
accession no. 1K4C; one subunit red) and an 
open mutant (Protein Data Bank accession no. 
3F5W; one subunit green) are shown, high- 
lighting the position of one likely proton sen- 
sor, H25 (not modeled in 3F5W). (B) MWC 
model with a single independent proton- 
binding site per KcsA subunit. Closed states 
(red) are indicated as "C" and open states 
(green) are indicated as "O," with subscripts 
indicating the number of protons bound to 
the channel. L a is the gating equilibrium in 
the absence of protons. K a doscd and Jf a open are 
the acid dissociation constants for the closed 
and open states, [tt ] is the proton concentra- 
tion. Each bound proton increases the gating 
equilibrium by a factor of y = X a dosed /iC° pen - 
(C) Two-site pH-sensor MWC model drawn 
as two tiers with protonation of site 1 along 
the x axis and protonation of site 2 along the 
y axis. Horizontal transitions depend on state- 
dependent pK a values as in B. Protonation of 
sites 1 and 2 increases the gating equilibrium 
(vertical transitions) by 7 = X al cl0K 7A al 0 P en or 
(J = ^ a2 closed /^a2 open , respectively. 
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thus shifting its closed-state pK a value to more acidic 
values. In contrast, a crystal structure of a mutant KcsA 
in a presumed open state (Fig. 4 A) showed a displace- 
ment of the TM2 helices, which would dramatically 
increase the solvent exposure of H25 (Cuello et al., 
2010b), yielding a higher pKa value, closer to 6.5. A simi- 
lar but less drastic change in pK a can be assumed for 
El 18, as it is located very close to H25 at the bundle 
crossing (~3 A in the closed crystal structure) and is 
less buried. 

The proposed model value for the wild-type intrinsic 
gating equilibrium suggested that KcsA is essentially 
fully closed in the absence of protons (L a = 2.5 x 10~ 12 ), 
which is also observed in electrophysiological record- 
ings (Cuello etal., 1998; Heginbotham et al., 1999). De- 
spite this remarkably stable closed state in the absence 
of protons, the modeled H25 sensor provides enough 
energy when fully protonated to significandy activate 
the channel by itself. If H25 is the pH sensor site 1 in 
Fig. 4 C, then full protonation of H25 raises the gating 
equilibrium constant to [CUol/tCio] with a value of 7 4 
Lo = 0.4, corresponding to a significant open probabil- 
ity, Po = 7 4 Lo/ {I + 7 4 Lo) = 0.29. The energy contrib- 
uted to gating by protonation of a single H25 residue is 
—3.8 kcal/mole in this model (Eqs. 4 and 5). Why then 
is the H25R channel completely closed rather than 
more open at high pH? The gating equilibrium L Q for 
the H25R channel is not 0.4 but is modeled as L Q = 7 x 
10~ 4 (Table 1). This discrepancy may result from an 
overestimation of the "strength" of H25 and neglect of 
unidentified pH-sensor residues that contribute to the 
activation of the channel. Alternatively, this difference 
may reflect the fact that arginine is an imperfect mimic 
for protonated histidine, and therefore the value for 
the gating equilibrium L a in the mutant channel can 
differ significandy from the value expected from the wild- 
type model parameters. Aside from the obvious steric 
differences between these residues, the charged guani- 
dinium group of an arginine residue may not localize to 
the same structural position as the charged imidazole 
group of histidine 25. For these reasons, we did not at- 
tempt to constrain our models such that pH-sensor mu- 
tations converted the energetic contributions of proton 
binding in wild type into a precisely corresponding 
change in L Q for the mutant channels. 

We modeled the secondary El 18 sensor as weak, with 
only a 1 .2-pH unit shift in pK a value between closed and 
open that fit the shallow proton dose-response of the 
H25R mutant (Fig. 1; Materials and methods). The en- 
ergy contributed to gating by protonation of a single 
El 18 residue is —1.6 kcal/mole in this model (Eqs. 4 
and 5). Therefore, El 18 cannot activate the wild-type 
channel on its own (L Q (3 4 = 1.5 x 10~ 7 ) . However, on the 
background of the H25R mutant that has an elevated 
L Q value compared with wild type, the weak secondary 
sensor can fully activate the channel (Fig. 1 B, red curve) . 



The pH-independent mutant H25R/E118A (Fig. 1 D) 
suggests that protonation of both H25 and El 18 are 
sufficient to account for the full activation of the pH- 
dependent gate in the control channel, assuming that 
the H25R/E118A precisely mimics the protonated H25 
and El 18. Alternatively, the mutations may have in- 
flated L Q because of steric or structural changes, resulting 
in an overestimation of the importance of these residues 
within the wild-type pH sensor. Notably, the assumption 
that the individual mutations H25R and E118A do not 
alter the remaining pH sensor but only alter L a did not 
accurately predict the L a value for the H25R/E118A 
double mutant (Eq. 6). In our model, £ 0 H25R/E118A = 3.1, 
corresponding to a Po value of M).75, was significantly 
different than the measured Po of ^0.98. Additionally, 
mutation of El 20 was assumed not to perturb the pK a 
values for El 18 or H25 to provide necessary constraints 
on the model parameters (see Materials and methods) . 
However, the model indicates that E120A increases L a 
when added to the H25R channel but slightly decreases 
L a when added to the H25R/E1 18A channel (Fig. 1 and 
Table 1 ) . These discrepancies suggest that the pK a val- 
ues might be somewhat altered by nearby mutations, 
and the total strength of each pH-sensor residue in the 
wild-type channel is therefore difficult to establish quan- 
titatively using the mutagenesis approach. In addition, 
other contributing residues have not been conclusively 
ruled out. 

What data supports the existence of other pH-sensor 
residues? Despite the dramatically pH-independent be- 
havior of H25R/E118A, other mutations of these sites 
have displayed pH-dependent activity. Three triple mu- 
tations of H25/E118/E120 that displayed pH-depen- 
dent activity were reported previously. First, H25R/ 
El 18A/E120A, a mutant that induced channel to chan- 
nel variability (Thompson et al., 2008) , displayed a small 
decrease in Po at higher pH values, indicating that the 
addition of E120A likely lowered the value of L Q in this 
triple mutant relative to H25R/E118A (Fig. 1 B). More 
convincingly, two other triple mutants maintained ro- 
bust pH-sensing capabilities. H25R/E1 18Q/E120Q was 
pH dependent with a Hill slope of ^2.5 (Thompson 
et al., 2008), indicating that the mutant had a similar 
sensitivity to protons as the H25R single mutation. Fi- 
nally, the H25A/E1 18A/E120A mutant was also pH de- 
pendent with a Hill slope of ~2.5. In both these latter 
mutants, the channel was completely closed at high pH 
and was fully activated by an increase in proton concen- 
tration, indicating the presence of an unidentified pH 
sensor (Fig. S2 and Table 1). 

If the KcsA pH sensor were distributed over a larger 
number of sites, the development of a realistic MWC 
kinetic model would be substantially more complicated 
as each additional site adds significantly to model com- 
plexity. For example, the case of three independent 
proton-binding sites increases the model to 250 states, 
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with seven free parameters describing the steady-state 
Po. On a conceptual level, however, a model includ- 
ing three or more independent proton-binding sites 
directly follows from the models presented, where 
each pH sensor has a state-dependent pK a value corre- 
sponding to its energetic contribution to gating the 
channel. Other sources of complexity are also possible, 
including interdependence between proton-binding 
sites, which we ignored in this study, as there was no ex- 
perimental result specifically indicating it. The physi- 
cal proximity of El 18 and H25 in the closed-state KcsA 
structure (~3 A), however, may suggest that proton 
binding at these residues (and thus the pK A values) 
may not be independent. 

Arginine residues modulate gating kinetics 
and influence proton binding 

In our model, mutation of residues near the pH sensor, 
such as the arginine residues, does not perturb the loca- 
tion of proton binding but must have a significant effect 
on either the pK a value or the gating equilibrium, L 0 , to 
alter pH-dependent gating. The pH dependence of the 
arginine mutants can be easily modeled by changing 
only L a (Table 1, model fit 1, and Fig. 2) or by changing 
both L Q and pK a values (Table 1, model fit 2, and Fig. 2) . 

When arginine residues were mutated in conjunction 
with E118A, an extremely dramatic change in gating 
kinetics occurred at intermediate pH values (Fig. 2 C) , 
suggesting an important interaction between El 18, 
R121, and R122 during gating. This interaction is also 
supported by the crystal structure of the El 18A mutant, 
where the positions of the arginines seem to be altered 
by the absence of the glutamate (Fig. 2 D). At interme- 
diate levels of activation, both double mutants displayed 
minute-long open periods interspersed with minute- 
long closings. This phenotype leads to an apparent 
channel to channel variability in pH-dependent open 
probabilities, although the pH dependence was clearly 
detectable (Fig. 2, A and B) . A macroscopic measure- 
ment, sensitive only to the ensemble average activity of 
many channels, would be blind to these effects but 
would likely record a shallower dose-response curve. 
These results suggest that the energy of an intermediate 
state is greatly increased relative to the energy of both 
the closed and open states in the two double mutants. 
In other words, the cooperative transition of subunits, 
fundamental to the MWC model, is altered or disrupted 
in these mutants relative to wild type. 

Conclusions 

In summary, we have shown that the pH dependence of 
KcsA arises in part from a histidine (H25) and a gluta- 
mate (El 18) located next to the closed bundle-crossing 
gate. A second glutamate (E120) does not appear to con- 
tribute to pH-dependent gating, and an apparent inter- 
action between E120 and a nearby arginine (R122) does 



not significantly contribute to pH sensing. We described 
the KcsA pH sensor using a two proton-binding site 
MWC model. The model offered a conceptual frame- 
work to rationalize the effects of mutagenesis and to 
propose thermodynamic contributions of proton bind- 
ing to channel gating. Histidine 25 is likely the central 
pH-sensor residue contributing the most energy to 
gating, whereas El 18 and possibly other unidentified 
residues probably enhance activation with smaller ener- 
getic contributions. R121 and R122 likely perturb the 
state-dependent pK a values for nearby pH-sensor resi- 
dues and also participate in setting the energy barrier 
crossed during the cooperative transition between closed 
and open. 
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